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Abstract 

Model-independent features of hybrid states with heavy quarks and 
adiabatic approximation for such systems are discussed. These general 
arguments are applied to the QCD-based constituent gluon model. 
The leading-order Hamiltonian is derived from the general formalism. 
The masses of heavy quarks hybrids, including charmonium hybrid 
states, are calculated. We predict ccg hybrids at 4.0±0.1 GeV and bbg 
ones at 10.5±0.1 GeV. Comparison with some other models are made. 
Model-independent features of hybrid states with heavy quarks and 
adiabatic approximation for such systems are discussed. These general 
arguments are applied to the QCD-based constituent gluon model. 
The leading-order Hamiltonian is derived from the general formalism. 
The masses of heavy quark hybrids, including charmonium hybrid 
states, are calculated. We predict ccg hybrids at 4.0±0.1 GeV and bbg 
ones at 10.5±0.1 GeV. Comparison with some other models are made. 
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The problem of exotic states with heavy quarks is widely discussed in 
the literature. Despite the fact that there is no unambigous evidence of ex- 
istence any exotics even for light quarks sector, intensive theoretical studies 
of heavy quarks exotics are carrying out (see [5] for review). The reason of 
this interest is that the heavy flavours spectroscopy is well understood, much 
more better than for light quarks sector. Presence of the slow nonrelativis- 
tic quark subsystem allows one to use the adiabatic (Born-Oppenheimer) 
approximation. 

So, there are slow quark subsystem and fast gluonic subsystem . The 
only mass parameter for gluonic subsystem is string tension a. Its energy in 
the order of magnitude is E g \ oc y/a; radius is r g i oc l/y/a; and velocity is 
f g i oc 1. Assume that quarks are confined by oscillator potential produced 
by gluonic subsystem . Then for quark subsystem we can estimate : 

E q g oc a 3/4 /y/E; r-gg oc 1/ W ma 3 / 2 (1) 
Hence , ratio of the periods of slow and fast subsystems equals : 

Tgg/T gl OC E gl /Egg OC (2) 

and radii ratio is : 

r q q/r g i oc A \j^a/m (3) 

It's easy to see that a typical small parameter is \J y/a/m. Unfortunately, 
because of the numerical reasons estimation ( 3 ) is not satisfied for c- and 
b-quarks . But the main point we stress is that the oscillator potential 
gives us quark subsystem energy of order yfo\j y/a/m and any unharmonic 
corrections will be of the next orders. 

So , we propose the following description of hybrids with heavy quarks . 
Suppose we have a Hamiltonian for qqg- system . We'll obtain the leading 
approximation in the limit r qq - — > 0. We can calculate gluon subsystem energy 
solving Schrodinger-type equation with our Hamiltonian . To obtain quark 
subsystem energy we ought to 1) substract from our Hamiltonian all terms 
proportional to the square of the interquark radius. 2) Averaged over gluon 
wave function these pieces of the full Hamiltonian produce us Hamiltonian of 
the quark subsystem which in turn gives us quark corrections to the energy 
of the gluonic subsystem . 
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Our starting point is the constituent gluon model described in [1] . Using 
Feynman-Schwinger representation for Green function one can write effec- 
tive Lagrangian for qqg system in the Euclidean space-time (spin degrees of 
freedom are neglected): 



2 fj, I I A 

(4) 

l\ = (fti + (1 - P^xf - \(f3i(zn) + (1 - f3i)(zn)) 2 

r 1 

l\ = {(3 2 z + (1 - (3 2 )x) 2 - \(fo(5r 2 ) + (1 - (3 2 )(zr 2 )) 2 

(We consider equal quark masses case, Coulomb interaction omitted). 
Here \i and \i g are effective masses of quark and gluon respectively ; z, z and 
x are quark, antiquark and gluon 3-coordinates ; T\ = z — x; r 2 = z — x. 

To deal with the square root terms we use the auxiliary field formalism 
[2]. Following the method of [2] we obtain the Lagrangian: 



2i i 

L = — + n + -^(1 + x 2 ) + /i(i 2 + + dp 1 iy 1 (p 1 )(l + Z 2 ) + 

jJL Zi JO 

+ -M/ 1 -^| T + (1^2))] (5) 
Jo vi{pi) 

To eliminate center-of -mass motion we can introduce variables (analo- 
gously to the [2]) : 

r li = z i -x i ; r 2i = z i -x i \ ^ = (1-6-6)^ + 6^ + 6^ (6) 

where £ 7 (7 = 1,2) are parameters to be defined from the condition that all 
terms of the type (Rfi) and (Rf 2 ) cancel ( but not (Rri) (rri) type terms!). 
It's easy to see that 

£ = Vi + SodPiVi(Pi)Pi t = (1 <_> 2 ) (7) 

— * 

After that the integration over center-of-mass coordinate R can be carried 
out . The resulting Lagrangian reads : 
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771 ^ A/f 

L = — + y + ^-[ /U (l-2a + 2e i 2 )+Mi 2 + ci+a 2 e i 2 ] + (l - 2)-Ma6(nr 2 ) + 
rfi(T 2 r? Ci (riri) 2 , , 1 , 1 rr9/ , 9 . . , (nr 2 ) 2 

+^ - y ^T" + ( ~ } " 2A { # ?(nri) (ai + M + a2 H^ } " 

-26M6 1 (a 1 + r))(hri)(f 2 ri) + £ 2 2 Mai(ai + r 1 ){r 2 r l ) 2 } + 

+ (1 <-> 2) + [Mfin) - 6«i(r 2 r 1 )][6 2 (f 2 r 2 ) - W^)] ^ ^ } (8) 

r x r 2 

where 

d/3 



Jo Jo Jo Jo 



\ r l x r 2\ 

i] = iii + /i 2 + fx g ; M = r] + ai + a 2 ; A = Mr] + aia 2 ^ — (9) 

r x r 2 

We assume that v ~ yfa. Expanding (8) on the y/a/m and using new 
variables 

Zi + Zi r u + r 2i _ 
n = — x 1 = ; pi = Zi - Zi = r u - r 2i (10) 

it's easy to see that in the leading order relative quark motion (with reduced 
mass ix 1 2 ) is decoupled : 

\^-r 2 ) 2 = \p 2 (11) 



In the leading order we have for gluonic subsystem : 
^ = 17(A f 3+ n i+ n 2)--^ : — \ 7T — ^- + ^{Hg+ai+a2+o- 2 rldi+a 2 rld 2 ) 

Z Z T'y Z V<2 Z 



(12) 
where 

n = c-6 + a/4= / d{3v({3)({3- l) 2 (13) 
•/ 

The kinetic term in the Lagrangian (12) is quadratic in velocities. Hence 
one can easily obtain Hamiltonian and all terms proportional to p 2 can be 
separated. The result reads : 
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f(0) = J_(Pjf L gl 

2fi g r 2 2(^ + n! + n 2 )r 2 

We assumed S-wave for quarks in ( 14 ) : 

< /HP S >= ^ (15) 

Averaging (14) over chosen gluon wave function we obtain a contribution 
to the oscillator potential for quark subsystem. In fact, the main contribution 
arise from the "potential" term of Lagrangian (12) and (14) is only a small 
correction. 

Zeroth approximation Hamiltonian of gluonic subsystem in Minkowski 
space is: 

H gl = + - Lg * + f 1 d(3u((3) + ^ + aV f 1 (16) 
9 2/i g 2(2U + fi g )r 2 Jo 1 v 11 2 Jo u(/3) V ; 

Here L g i is the angular momentum of the gluon ; p r is the radial momentum. 
We believe that physical states of gluon are electric and magnetic ( see [1]). 

^2 

For lowest states < L gl >= 2. 

Here we are to use some variational procedure. We choose an anzats for 
the effective string energy v{(3) and define free parameters from the mini- 
mization of gluon subsystem energy . After that quark correction with this 
ansatz is calculated from (11), (12) and (14) . The simplest ansatzes for v(f3) 
are : 

u((3) = Const (Ansatzl) 
u(p) = d + C 2 (f3 - l) 2 (AnsatzII) 
u((3) = or (Ansatzl II) 

Without Coulomb interaction it's worthwhile to write out values of inter- 
est in the dimensionless form : 



Egi/yfa = 4.30 ± 0.05; E qq /{y/a^ y/a/m) = 1.4 ± 0.2 (17) 
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These values were obtained by means of trial radial wave function of the 
Gaussian type. Numerically gluon energy is about 1.7 GeV , quark correction 
value equals hundreds of MeV. 

Three Coulomb terms should be added to the full Hamiltonian: 

V coul = -—^(- + -) + -- 18 

2 ri r 2 6 p 

Since a s /6 <C 1 quark- ant iquark repulsion gives only small contribution into 
the resulting energy . 

But the main problem as it is usual for most of the (potential-like) models 
is to define absolute scale of mass. Some additional constants are needed to 
fit hadronic spestrum and these constants vary for different sectors of the 
theory. Following the recipe described in [ 1 ] we use the constant for heavy- 
light sector multiplied by two ( cause we have two strings vs. one string for 
heavy-light meson , and one can expect that perimeter-type terms for Wilson 
loop or hadronic shifts could be twice larger than for the only string). The 
procedure used is : 1) to define constant for heavy-light meson for given heavy 
quark mass (light quark is massless) . By means of variational calculations 
with known masses of the heavy-light S-states (1975 Mev for D-meson and 
5300 MeV for B-meson - accuracy more than tens of MeV is not needed) we 
define additive negative constant Co- 2)The resulting mass is: 

Mggg = 2m , + E gl + Egg + 2C Q ( 19) 

The main results of this paper are summarized in Table I . There hybrid 
masses for various heavy quark masses are listed. The parameters are a s = 
0.3; a = 0.18GeV 2 . In the last column results obtained by method of [ 1 ] 
are shown. 
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Table I. 

Calculated hybrid masses, GeV 



m r = 





Ansatz I 


Ansatz II 


Ansatz III 


Potential 


1.2 GeV 


4.008 


3.951 


3.873 


4.089 


1.3 


4.004 


3.95 


3.875 


4.082 


1.5 


3.997 


3.946 


3.878 


4.069 


1.7 


3.991 


3.943 


3.88 


4.058 


4.2 


10.828 


10.810 


10.647 


10.647 


4.5 


10.551 


10.521 


10.491 


10.585 


4.9 


10.55 


10.521 


10.492 


10.583 


5.2 


10.548 


10.519 


10.491 


10.58 



m b -- 



From Table I one can conclude that our model anticipate for c.o.m. of 
hybrid multiplet 

J 4.0±0.1GeV for ccg 

\ 10.5±0.1GeV for bbg { ' 

Quantum numbers of hybrids are given in [ 1 ] . 

Spin-dependent forces are not considered in this paper ( the most impor- 
tant are (LS). They will be discussed in the forthcoming studies. 

The result (19) is in a good agreement with that ones of [3] and [4]. 
Unfortunately, this is a rather accident fact cause in [3] , [4] and in this paper 
various prescriptions for additive constants in the interquark potential were 
used. Thus we believe that there is no reasons for exceed optimism and 
additional attempts should be made to understand the origin of the constant 
term. 
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